Experimental Polarization State Tomography using Optimal Polarimeters 
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We report on the experimental implementation of a polarimeter based on a scheme known to 
be optimal for obtaining the polarization vector of ensembles of spin-| quantum systems, and the 
alignment procedure for this polarimeter. We also show how to use this polarimeter to estimate the 
polarization state for identically prepared ensembles of single photons and photon pairs and extend 
the method to obtain the density matrix for generic multi-photon states. State reconstruction 
and performance of the polarimeter is illustrated by actual measurements on identically prepared 
ensembles of single photons and polarization entangled photon pairs. 



I. INTRODUCTION 

Many promising applications in Quantum Information 
(e.g. Quantum Computation and Quantum Communi- 
cation) demand accurate state estimation. For many of 
them it is compelling to implement state estimation tech- 
niques that are both fast and consume as few copies of the 
state as possible. Research in improving the efficiency of 
quantum state estimation techniques is an area of active 
theoretical study 0, IE SS IE 01 with much focus on 
spin-i systems (qubits). Experimental reports on state 
estimation are fewer [Sj, partly because many schemes 
call for a joint measurement on an ensemble of qubits 
which is not always possible to implement. 

One way to realise qubits experimentally is to use 
single photons and consider their polarization degree of 
freedom, as this can be described by a two dimensional 
Hilbert space. Estimating the polarization state of an en- 
semble of single photons (called polarimetry in classical 
optics) is equivalent to estimating the state of the qubit 
ensemble. This makes efficient polarimetry interesting to 
Quantum Information. 

Polarimetry that uses the least number of measure- 
ment outcomes is said to be minimal. Minimal polarime- 
try techniques in classical optics have been known for 
a long time and a lot of work in their optimization has 
been done |E EE 0> 1^ • While these classical meth- 
ods perform well in estimating the polarization state for 
single photon ensembles in the limit of large numbers, 
their performance in the regime of extremely low light 
intensity (single photon level) was uncertain and it was 
not obvious how to use them in estimating non-classical 
states of light. For this reason, progress in polarimetry 
at the single photon limit will assist in many areas, in- 
cluding characterization of faint sources of light, classical 
ellipsomet ry l9 || , advanced quantum key distribution pro- 
tocols 10, Il5lll6| and studies of the fundamental aspects 
of quantum theory ^3 ■ 

In discriminating between different estimation tech- 
niques we distinguish between methods that are minimal 
and those that are minimal and optimal [E ■ optimal 
methods have the best asymptotic efficiency in determin- 
ing an unknown state when averaged over all possible in- 



put states. This can be used for an operational definition 
of minimal and optimal state estimation for ensembles of 
prepared quantum systems. It is the technique that pro- 
vides the best improvement to our estimated state for 
each additonal copy taken from the ensemble. Recently 
Rehacek et al. proposed such a method for state estima- 
tion of polarization based single qubits , which can be 
viewed as an extension of classical techniques 0, 0| 
to quantum systems. 

In this paper we address the experimental problem of 
implementing the optimal state estimation method de- 
scribed in [7| by using a complete four output polarime- 
ter with no moving parts. We describe that polarimeter 
in section II by reviewing the theory of optimal polariza- 
tion state estimation and explain our implementation. In 
section III polarization state estimation of multi-photon 
states is addressed. In section IV we elaborate on the 
alignment procedure to make the polarimeter perform 
optimally. Experimental state reconstruction on ensem- 
bles of single photons and photon pairs with high fidelity 
will be illustrated in sections V and VI. 



II. STATE ESTIMATION USING THE OPTIMAL 
POLARIMETER 

The polarisation state of light can be characterized 
using the three Stokes parameters 81,82,8^, possibly 
augmented by an intensity 8m- Together they form a 
Stokes vector 8 and, when normalized, it is written as 
8 = {1, 81/ 8rm 82/ 8m, 83/ 8m)- A rcditced Stokes vector 
8r ~ {81, 82, 83) /8m identifies a point in the Poincare 
sphere (in this paper reduced Stokes vectors are denoted 
by an r subscript) [l^ . 

A minimal scheme of estimating the Stokes vector re- 
quires exactly four detector readings, which corresponds 
to finding the overlap of the unknown Stokes vector with 
four non-coplanar vectors that define a tetrahedron in the 
Poincare sphere (Figure^]) . These four non-coplanar vec- 
tors define four measurement operators Bj that govern 
the detector readings and form a set of complete Pos- 
itive Operator Value Measurements (POVM) The 
tetrahedron geometry defines the largest volume that can 
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HWP at 22.5° 



FIG. 1: Four (reduced) Stokes vectors in the Poincare sphere 
that form a tetrahedron define the optimal POVM operators 
used for polarization state estimation. The tetrahedron gives 
the largest volume encompassable by a vector quartet in the 
sphere making it the optimal measurement when using four 
POVMs. 



be enclosed by a vector quartet in the Poincare sphere, 
making it the optimal estimation technique when using 
four POVMs [lll|2l|. Such a state estimation technique 
is also unbiased in the asymptotic hmit because the total 
distance of any vector in the Poincare sphere to all four 
POVM vectors depends only on the vector's magnitude. 
In other words, the orientation of the unknown vector 
does not affect the accuracy with which it is estimated 



We shall denote the tetrahedron's reduced Stokes vec- 
tors by &ir, ^2rj ^3ri ^4r shown in Figure n and write 
their corresponding normalized vectors as 6i, 621^31 
Each measurement operator Bj may be expressed as 



(1) 



where a — (cto, fi , (72 , ca ) , (Tq being the unit matrix and 
fi,2,3 the Pauli matrices. 

In an experiment we associate each operator Bj with 
a detector bj. The average intensity falling on detec- 
tor bj is denoted as Ij . Thus expectation values of the 
tetrahedron operators are related to detected intensities 
as follows: 



^ = (B., 
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with 



It^Y.^^ (2) 



Writing the intensities as a vector I = {Ii, I2, 13, h)/ It 
gives us the Stokes vector 



n-s ^ s = 11- 



(3) 



where 11 is referred to as the instrument matrix. Each 
row of this matrix is composed from a vector bj. One 
possible instrument matrix of the ideal polarimeter is: 
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Light of 
unknown („) 
polarization ^ 



PPBS, \yp/ \ PBS 



i_ Quartz 
Plates 



QWP at 45° 
b4a^PBS 

FIG. 2: Practical implementation of the tetrahedron po- 
larimeter that achieves the ideal instrument matrix. Each 
detector hj is associated to the tetrahedron vector hj. The 
partially polarizing beam splitter (PPBS) separates incom- 
ing light according to polarization, and quartz plates remove 
unwanted phase shifts. Light leaving the PPBS is passed 
through waveplates and polarizing beam splitters (PBS) to 
be projected on two different bases (±45° basis for transmit- 
ted light and the circular basis for reflected light) . 



Experimental realisation of this instrument matrix is 
achieved by the polarimeter shown in Figure |2 The first 
component of the polarimeter is a partially polarising 
beam splitter (PPBS) that has a particular amplitude 
splitting ratio for incoming light, most easily determined 
using Jones vector notation for polarization. The am- 
plitude division coefficients of the PPBS x and y obey 
energy conservation -|- = 1. The PPBS takes 



horizontally polarized (H) light (p) to the polarizations 
(g) and (q) in the transmitted and reflected arm, re- 
spectively, and vertically (V) polarized light (") to (°) in 

transmission and (^) in reflection. 

If we project light in the transmitted arm of the PPBS 
on the ±45° polarization basis and light in the reflected 
arm onto the Left/Right circular polarization basis, the 
tetrahedral arrangement of the vectors hj are ensured 
with the following relations: 



2V3 



and y 



(5) 



Detailed steps are given in Appendix A. 

Partially polarized light can be described using a den- 
sity matrix (or coherency matrix). If we write the entries 
of the density matrix as a column vector p, we can deter- 
mine them from the Stokes vector 22] using the following 
transformation: 

P = \^i-s 



/I 1 \ 
1 i 
1 -i 

V 1 -1 / 



• s 



(6) 



The columns of the matrix Ti are the Pauli operators 
written as column vectors Fi = [ctq, di, 02, 03). The ma- 
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trices B ^ and Fi can be combined into a single matrix 
T:=irin-i ^ p^T-I (7) 

which might be referred to as a tomography matrix as 
it directly relates the detected intensities to the density 
matrix of the state. 

III. POLARIZATION STATE TOMOGRAPHY 
FOR ENSEMBLES WITH MULTI-PHOTONS 

The instrument matrix scheme above can be extended 
to perform polarization state tomography on ensembles 
of multi-photon states. James et al. |8j have described 
a similar state estimation method. We follow their ap- 
proach but use our optimal and instrumentally motivated 
measurement operators, thereby reducing any ambiguity 
over the choice of operators. 

The simplest multi-photon system is a photon pair de- 
tected by testing for coincidence in the detection time of 
their component photons. In our measurement process 
each member of the photon pair is passed through a po- 
larimeter. Given two polarimeters 1 and 2, each with four 
detectors b^^ and b^j, respectively, (ii,i2 = 0, 1,2,3), we 
will have 16 possible coincidence combinations between 
the detectors (see Figure |2J . Each coincidence rate is 
governed by an operator composed from the individual 
detectors' measurement operators. If we denote again the 
measurement operator of detectors bi^ and hi^ as Bi-^ and 
Bi^ , and the coincidence count between them as Ci-^^i^ , we 
can express the coincidence rates as a linear function of 
a 2-photon polarization state vector S2'- 

^ = = {lbl<g>hl)-S2, (8) 

ct 4 4 

4 

with Ct = ^ Cii,i2 

Here, 5*2 is the Stokes vector equivalent for a 2-photon 
system and Cj is the total number of observed coin- 
cidences. We now have the set of measurement opera- 
tors governing the coincidence pattern. The sixteen co- 
incidences Cij.ij can be written in column vector format 
C2 =(ci,i, C1.2, 04^4). If we define the 2-polarimeter in- 
strument matrix as 112, we obtain an instrument response 
analogous to (3): 

C2 = n2 • ^^2 ^ ^2 = 112' ■ C2 (9) 

Thus we obtain the density matrix of the 2-photon 
state by constructing the analogous 2-photon expression 
for equation ®: 

P2 = • ^2 = ^2 • C2 (10) 



detector outputs 





photon 1 


"11 


Multiple 

coincidence 

unit 


Source 
of N 

correlated 
photons 


PI 


b2i 




b3i 


b4i 


photon 2 


bi2 


P2 


b22 


photon A/ 


b32 


b42 




bl N 


PW 


b2A/ 




b3N 


b4A/ 



FIG. 3: Scheme for estimating the polarization state of 
an ensemble of A^-correlated photons using A'' polarimeters 
(P1,P2,...,PA'^). A multiple coincidence circuit identifies the 
4^ possible coincidence combinations. For photon pairs 
{N—2), two polarimeters are used giving 16 possible coinci- 
dence combinations. Several copies of the state are processed 
giving a coincidence pattern used in estimating the polariza- 
tion state of the ensemble. 

Each column of r2 is the product of two Pauli operators 
Cii (8) CTia *2 — 0, 1, 2, 3) written in column vector for- 
mat and T2 is the tomography matrix for the 2-photon 
state. 

It is now straightforward to generalize this concept to 
obtain the density matrix for states of N correlated pho- 
tons. Using N polarimeters, we obtain the pattern of 
TV-fold coincidences to build up the coincidence vector 
Cn which is used to find the iV-photon Stokes vector 
and density matrix: 

Sn = • Cn, (11) 
Pn — '^^N ■ Sn — Tn ■ Cn (12) 

Each row of the instrument matrix IIjv is given by 
( 3 ''ii ® i ■ • • ® i ^iw ) each column of Fjv is the prod- 
uct of N Pauli matrices ^ ai^... ^ {in = 0, 1, 2, 3 
and n — 1,2, N). This generalized approach will 
work for all four-detector polarimeters in multi-photon 
analysis schemes (Figure I^J. We note that — although 
the polarimeters are optimal for estimating single qubit 
states — it is an open question if the scheme above is 
optimal in estimating multi-photon systems. 

IV. PHASE CORRECTION AND POLARIMETER 
CALIBRATION 

IV a. Removing Unwanted Phase Shifts 

In the presented polarimeter, an ideal PPBS has the 
nominal beam splitting ratio © and also rotates the po- 
larization state of light leaving the beamsplitter into the 
correct polarization basis ^^l- Such beamsplitters, how- 
ever, are not easily available and their design is the focus 



4 




HWP Angle (deg) HWP Angle (deg) HWP Angle (deg) 



FIG. 4: Instrument response of the polarimeter to linearly polarized light. The data points show the variation in relative 
intensity at each detector with respect to the angle of the half wave plate (HWP) in the polarisation state preparation. The 
solid lines show the expected intensity modulation for an ideal device for each HWP setting 1131 . scaled for appropriate 
detector efficiencies. Error bars are smaller than the point markers. Panel (a) shows the relative intensity at detector 2 without 
compensation plates. Panels (b) and (c) are taken with compensation for phase shifts. The oscillation in (a) is out of phase 
and also of lower amplitude compared to the phase corrected behaviour of detector 2 in panel (b). Plots in (c) shows a lower 
amplitude because light in the reflected arm is not projected on a linear polarization basis. 



of active research |23|. We therefore use beamsplitters 
with only the nominal intensity splitting ratio. 

A PPBS without phase shift diverts light in state —bjr 
(that is conjugate to a tetrahedron vector bjr) from detec- 
tor bj. General beamsplitters, however, lack this phase 
preserving property. The result is that input of conju- 
gate states —bjr does not stop light from reaching the 
associated detectors. This suggests an easy alignment 
method for correcting any unwanted phase shifts with 
birefringent compensation plates. 

For phase correction we first prepare high quality H- 
polarized light using polarizers of extinction ratio 10^. 
With one subsequent half wave plate (HWP) and one 
quarter wave plate (QWP) we can then prepare any po- 
larization state on the surface of the Poincare sphere. 
Compensator plates (0.5mm thick quartz) mounted on 
rotating stages were placed at each output arm of the 
PPBS, and light with a conjugate polarization state was 
sent to the polarimeter. For each polarization state —bjr 
the compensator in the relevant output arm was rotated 
until the detector hj received no light. Two input states 
(one for each output arm) were sufficient to compensate 
for the unwanted phase shifts. 

We verified the compensated polarimeter behavior 
with linearly polarized light prepared using only the po- 
larizer and HWP (this reduces preparation errors due to 
residual errors in the QWP). The prepared states have 
a Stokes vector of the form (l,cos4'0,sin4V',O), where ip 
is the angle of the HWP, so the normalized response of 
detector 1, for example, will be 

/i = l + -y/icos4V'-l-y^sin4V'. (13) 

Passively quenched silicon avalanche photodiodes were 
used as detectors allowing us to perform photon counting. 
The number of photons accumulated at each detector 



output was noted for each angle of the HWP. The results 
are shown in Figure^ 

The results show that the response of the compen- 
sated polarimeter is very close to ideal. The extrema 
of our measured intensities are less than 1° (of HWP 
angle) away from their nominal positions. This means 
that our actual measurement vectors are pointing in the 
same direction as the ideal tetrahedron vectors, although 
their magnitudes will be different due to imbalanced de- 
tection efficiencies. While this renders the asymptotic 
efficiency of our polarimeter less than ideal, it still rep- 
resents the optimal setup for the collection efficiencies 
we can achieve. In other words, we are maximizing the 
volume defined by our experimental POVM vectors [2H . 

This measurement result is limited by the accuracy of 
our rotation controllers. Our waveplates are mounted on 
rotary motors with an accuracy of 0.3°. The polarizing 
beam splitters in the output arms have an extinction ra- 
tio of 10^ and the waveplates' optical path length differ 
from their nominal values by less than 2%. 



IV b. Calibrating the Polarimeter 

We calibrate the instrument matrix of this polarimeter 
to account for all residual phase shifts and coupling ineffi- 
ciences. A general calibration technique for four detector 
polarimeters ( "equator-poles method" ) was described by 
Azzam et al. |l8l |. Incidentally, the phase dependency 
measurement shown in Figure 01 was an essential part of 
this calibration. 

Using this technique we are able to find the correction 
terms needed to be made to our ideal instrument matrix. 
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Fidelity of Estimated States 




90 1 80 270 
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FIG. 5: A set of polarization states {Si) equally distributed 
over the Poincare sphere surface was generated; photons from 
each of these states were sent to the polarimeter, from which 
an estimated state (Se) is obtained. The figure shows the 
fidelity of the estimated state to the prepared state, ^(Se-Si). 
It is roughly constant over the Poincare sphere, showing that 
the polarimeter is an unbiased polarization state estimator. 



A typical corrected instrument matrix is shown below: 



He 



0.962 1.051^/1 



0.991 1.031 



1.010 
1.032 



-1.045ji 



-1.009 



0.920 
-0.956 
0.005 
0.029 



0.005 \ 
-0.005 



-O.945JI 



1.011 



The uncertainty for each of the correction terms above 
is on the order of 0.002. We see that the deviation from 
entries in the ideal instrument matrix is on the order 
of a few percent. 

The phase correction and calibration steps presented 
above must take into account the wavelength of the in- 
put light because optical elements are specified to per- 
form only within a certain bandwidth. The polarimeter 
was built to study the polarization state of light coming 
from a spontaneous parametric down conversion (SPDC) 
source [2^ with a spectral bandwidth of 4.740±0.014 nm 
centred around 702nm. The same light source was used 
for phase correction and polarimeter calibration and the 
experiments described in the remaining sections. 



V. EXPERIMENTAL STATE TOMOGRAPHY 
FOR ENSEMBLES OF SINGLE PHOTONS 

The ability of the tetrahedron polarimeter to estimate 
polarization states without bias was tested by preparing 
a set of pure polarization states equally distributed over 
the Poincare sphere. In this way we could better identify 
regions that suffer poor state estimation (if any). 

Computer controlled motors were used to rotate wave- 
plates (after a H-filter) in preparing the set of polariza- 
tion states. The Stokes vector of a pure polarization state 
can be expressed as 5 = ( 1 , cos 26 cos 2[ip + S], — cos 2[ip + 



5] sin 25, — sin 2 -I- 5] ) , where S and ip are the QWP and 
HWP angles, respectively. Thus any set of coordinates 
(characterized by the polar and azimuthal angles) on the 
Poincare sphere can be expressed in terms of the wave- 
plate angles. 

For each set of angles, the detectors accumulated pho- 
ton detection events for one second giving a particular 
vector / from which an estimated Stokes vector Se and 
probability density matrix pe can be obtained via equa- 
tions (PJ and To calculate the distance of the esti- 
mated state from the (ideal) prepared state pi {Si), we 
use the Uhlmann fidelity, defined as (tr[y\/pip7\/Pi])^ 
j25[ E^ . For pure states this quantity reduces to the 
overlap of their Stokes vectors ■^{Si ■ Sg). 

The fidelity was mapped to the appropriate polar and 
azimuthal coordinates on the Poincare sphere (FigureEJ; 
linear polarization states correspond to a polar angle of 
0°. The average fidelity for the whole map is 99.8% with a 
minimum fidelity of 98.4±0.9% (the cumulative photon 
count per point is approximately 2000). There are no 
systematic areas of low fidelity even when wedge errors 
in the state preparation waveplates cause count rates to 
drop. This indicates that the polarimeter estimates all 
pure polarization states equally well. 

Fidelity does not distinguish between errors introduced 
in state preparation from errors in the state estimation 
process. Therefore we have characterized our state prepa- 
ration apparatus independently and are confident that 
their contribution to the error in calculated fidelity above 
is on the order of ±0.01%. Thus we assign the residual 
difference in fidelity to imperfections in the detection ap- 
paratus. 



VI. EXPERIMENTAL STATE TOMOGRAPHY 
FOR A TWO PHOTON ENSEMBLE 

We will now illustrate how to use two polarimeters 
to perform polarization state tomography on a 2-photon 
state generated from an SPDC source. First, two po- 
larimeters were correctly aligned and after calibration 
their instrument matrices were found to be: 

0.9997^/1 -0.041 \ 

0.0602 



and 



0.903 0.927 



1.124 1.135ji -1.014/2 



0.995 
0.978 



-1.079 
-0.983 



0.001 
0.003 



0.913, /§ 



-0.936 



1.074 1.171^/1 



0.983 0.8804 
1.082 -1.172 
0.862 -0.^ 



0.913 
-1.044 
0.001 
-0.002 



-0.082 \ 
0.004 
-0.9625^1 



0.867, 
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(a) Real Part (b) Imaginary Part 




FIG. 6: The density matrix of a Bell state obtained 
by linear reconstruction from photon pairs. 

We then arranged for the SPDC source to generate 
photon pairs that are detected as a maximally entan- 
gled Bell state 1^'"'"). Bell states created via SPDC are 
typically characterized by a polarization correlation ex- 
periment, from which a visibility value can be obtained 
|2J|. The visibility measured in the HV and ±45° basis 
was above 97.7 ± 2%; such a high value is usually taken 
as evidence of a high degree of entanglement. 

The photon pairs were passed through the polarime- 
ters and the pattern of coincidences between them was 
observed. The 16 observed coincidence rates (collected 
using the scheme similar to [27|) make up the coinci- 
dence vector G =(21444, 1505, 24104, 26002, 979, 24716, 
23210, 22447, 21661, 30752, 24061, 268, 19010, 23692, 
339, 17695). 

Using this vector with equations and (|10|l we obtain 
the density matrix whose real components are 

/ -0.002 -0.01 -0.03 -0.024 \ 
. , -0.01 0.506 0.485 0.025 

\ -0.024 -0.024 0.009 -0.003 / 

while the magnitude of the imaginary components are 
below a value of 0.04 (see Figure 0. 

The uncertainty in each of the above terms is on the 
order of 0.011. The Uhlmann fidelity of this state to the 
ideal state was found to be 0.990 ± 0.014. Error 

bars in all cases were computed by numerical derivation 
and propagated Poissonian counting noise. The propa- 
gated error bars result in an estimated density matrix 
compatible with the ideal |^'+)(^'+| state. 

VII. CONCLUSION 

In this paper we have illustrated a simple align- 
ment procedure for optimizing the tetrahedron polarime- 
ter. Phase shifts introduced by a commercially avail- 
able PPBS were easily corrected by phase compensation 
plates in the output arms. The response of the compen- 
sated polarimeter was measured over a dense sampling 



of states on the Poincare sphere, and found to be similar 
to that of an ideal device. This shows that beamsplit- 
ters need only have the nominal intensity splitting ratio, 
making optimal polarimeters more accessible. 

We also described an instrument ally motivated method 
for constructing the measurement operators governing 
light distribution to each output of the polarimeter. This 
instrument based approach also allows a convenient gen- 
eralization to obtain measurement operators governing 
multi-photon coincidences. These operators can then 
be applied to the linear reconstruction of multi-photon 
Stokes vectors and their density matrices. 

Optimal polarimeters were then used for estimating 
the polarization state of experimentally prepared ensem- 
bles of single photons and photon pairs in a Bell state. 
The estimated states were evaluated by computing their 
fidelity to the (ideal) prepared states. We found an aver- 
age fidelity above 99.8% in all our experiments. Thus we 
have built and demonstrated the use of optimal 4-output 
polarimeters in multi-photon polarization state tomogra- 
phy. 

While preparing this document, it came to our atten- 
tion that a similar four output polarimeter was suggested 
independently in a recent paper . 
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APPENDIX A 

To derive the intensity splitting ratio of the PPBS, we 
first express the polarization states are expressed using 
Jones vectors unless we are describing the tetrahedron 
vectors hj. The tetrahedron (Stokes) vectors have 
the scalar product property 

h-'^k = \ + \5jk- (14) 

Recalling the parameters of the intensity splitting ra- 
tio of the PPBS X and y we see that a general input 
polarization state (^) leads to the polarizations (^^) and 
(^^) in the transmitted and reflected arms of the PPBS 
respectively. In our polarimeter, light leaving the arms of 
the PPBS must be analyzed in two different polarization 
bases. Two orthogonal vectors that form a basis may be 

expressed as and ("""'^J'"^)- This leads for ex- 

ample to the normalized light intensity falling on detector 
bi 

/i//t = |aa;cos6'-|-/3^/e-"^sin6'|^ (15) 
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We choose a different measurement basis for detectors 
3 and 4, for example light reaching detector 3 is 



ay cos ( 



j3xe 



-'■^'sine*' 



(16) 



Using the vector hi as an example, equation (0) allows 
us to express the operator Bi in terms of the measure- 
ment basis to fulfill equation H15|) : 



(Bi) 



\ax cos ( 



(iye 



smf 



X cos a ye 



(17) 



(18) 



The following choice of Bi fulfills this condition 

xcosO 
y%j sin 9e^ 



Bi = 



X cos f 



ysmt 



) ■ (19) 



Since the tetrahedron can be oriented arbitrarily we 
choose for convenience to measure the 45° linear polar- 
ization basis {0 = tt/A, (/) = 0) in the transmitted arm 
and the circular polarization basis {6' = tt/4 0' = 7r/2) 
in the reflected arm. This reduces the measurement op- 
erators to only the beamsplitting parameters x and y 



Bi 1 1 ±xy 
B2 ) 2 \ ±2;?; 



B3 



1 ^ ^^y 

2 \ zLixy x^ 



which together with equation Q allows us to express all 
tetrahedron vectors in terms of x and y 




y — X 


V T2xy J 



From equation (|14() . we can write: 

^ ^ 2 - - 2 
61 • &2 = g and &i • 63 = -. 

This allows us to obtain an equation in x alone 
36a;® - 24a;* -I- 1 = 0. 



(20) 



(21) 



(22) 



The last equation gives two solution sets; we choose 
the set where a;^ = ^ H — ^= ^ y^ = h h=- 

2 2\/3 " 2 2^3 
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